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Quenches in an ion chain can create coherent superpositions of motional states across the linear- 
zigzag structural transition. The procedure has been described in [Phys. Rev. A 84, 063821 (2011)] 
and makes use of spin-dependent forces, so that a coherent superposition of the electronic states 
of one ion evolves into an entangled state between the chain's internal and external degrees of 
freedom. The properties of the crystalline state so generated are theoretically studied by means of 
Ramsey interferometry on one ion of the chain. An analytical expression for the visibility of the 
interferometric measurement is obtained for a chain of arbitrary number of ions and as a function 
of the time elapsed after the quench. Sufficiently close to the linear-zigzag instability the visibility 
decays very fast, but exhibits revivals at the period of oscillation of the mode that drives the 
structural instability. These revivals have a periodicity that is independent of the crystal size, and 
they signal the creation of entanglement by the quantum quench. 
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I. INTRODUCTION 

Quenches of quantum many-body systems provide im- 
portant information on the thermodynamic properties of 
physical objects close to phase transitions. They give 
insight into the statistical mechanics of closed systems, 
and can find applications for quantum information 
3 . Among several proposals discussed in the literature, 
some set-ups make use of the coupling with a quantum 
system, a spin, to drive a quantum phase transition in a 
larger physical object acting as environment. These dy- 
namics are associated with decay of spin coherence that 
has been shown to exhibit universal features [3H7]. 

In a recent article J5| these dynamics were studied in 
a system of trapped ions, when these form a linear ar- 
ray in a linear Paul trap [9]. Here, the spin is an internal 
transition of one ion of the chain, while the vibrational 
excitations of the chain itself acts as an environment. 
Realizing a Ramsey-type of interferometer with the in- 
ternal transition of the ion, a quench is performed by 
the mechanical effect of light associated with the absorp- 
tion and the emission of the laser photon. This quench 
is carried out when the chain is close to the mechanical 
instability at which it undergoes a transition to a zigzag 
structure [10] . It was shown that the visibility of the 
Ramsey interferometer, which is found by measuring the 
population of one internal state of the ion transition af- 
ter the second Ramsey pulse, allows one to access the 
autocorrelation function of the chain at criticality [5]. 

In Ref. |llj it was proposed to use the spin excita- 
tion to create a superposition of two different crystalline 
structures across the linear-zigzag structural transition. 
The superposition can be accessed by driving the elec- 
tronic transition of one ion of the chain in a set-up 
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FIG. 1. (Color online) A quench across the linear-zigzag in- 
stability is performed by exciting the central ion with a laser 
pulse in presence of spin-dependent forces. In (a) the collec- 
tive motion is initially in the ground state of a zigzag struc- 
ture, and the central ion in the internal state \g). A 7r/2 laser 
pulse prepares it in the superposition (\g) + |e))/v2. (b) A 
tighter state-dependent potential, acting only when the ion 
is in state |e), induces conditional dynamics such that the 
ions' internal and external degrees of freedom get entangled. 
When the state-dependent potential is sufficiently tight, the 
excited component will start oscillating around the equilib- 
rium positions of the linear chain, (c) A laser pulse performs 
a — 7r/2 rotation on the central ion. The final occupation 
of the ground state \g) as a function of the time between the 
two pulses contains information on the properties of the chain 
across the linear-zigzag instability. 

where an external field makes the trap frequency spin- 
dependent [TTJ[T2]. I n these settings, a first laser pulse 
prepares the ion in a coherent superposition of the elec- 
tronic states, which evolves into an entangled state be- 
tween the chain's internal and external degrees of free- 
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dom as sketched in Fig. [T] The properties of the crys- 
talline state so generated were studied by evaluating nu- 
merically the visibility after a second laser pulse is ap- 
plied, as shown in Fig. IJc). The visibility of the inter- 
fcrometric signal was shown to exhibit a fast decay, in 
agreement with studies performed in other settings [5l— 
[7], while for longer times quasi-periodic revivals of the 
visibility were found. 

In this paper we analyse the dependence of the signal 
visibility on the system parameters for the set-up pro- 
posed in Ref. [IT] . We determine analytically the expres- 
sion of the visibility and study its behaviour close to and 
across the classical linear-zigzag instability, for different 
numbers of ions. We find that the revivals observed in 
the visibility as a function of the time t elapsed after the 
quench are characterized by the frequency of the zigzag 
mode, and persist when the number of ions is increased. 
The analysis of the spectrum of the visibility signal as a 
function of t shows the presence of squeezing and entan- 
glement that are generated by the quantum quench. 

The article is organized as follows: In Sec. [IT] the pro- 
posal of Ref. [UJ is summarized. The theoretical model 
is presented in Sec. |III[ which also includes the detailed 
evaluation of the visibility signal. The behaviour of the 
visibility is analysed in Sec. |IV| and the conclusions are 
drawn in Sec. [V] Theoretical details for the derivation of 
the results in Sec. |III| are given in the appendices. 



II. 



RAMSEY INTERFEROMETRY WITH AN 
ION COULOMB CRYSTAL 



In this section we briefly review the physical model at 
the basis of this work. A string of N ions of mass m and 
charge q is confined in a trap, forming a zigzag structure 
close to the linear-zigzag mechanical instability. A two- 
level transition of the central ion is driven by two laser 
pulses separated by a time interval t and which perform 
a 7r/2 and — ir/2 rotation of the dipole, respectively. The 
pulses are short so that the crystal dynamics can be ne- 
glected during their duration [8] [13] . Under the assump- 
tion that both internal states of the dipolar transition 
are stable, a spin-dependent force is applied, such that 
the stable configuration of a finite chain is a linear struc- 
ture when the ion is in the excited state [TTJ. Therefore, 
during the time elapsed between the two pulses, the crys- 
talline state undergoes conditional dynamics dependent 
on the internal state, that lead to entanglement between 
internal and external degrees of freedom [8] [IT] . 

In the following we denote by \g) and |e) the two in- 
ternal states of the central ion, and omit to write the 
internal state of the other N — 1 ions since this remains 
unchanged. Before the first pulse, the state of the central 
ion and crystal motion reads \ip{fy) = \g) |</>(0)) where 
\4>(0)) can be either the ground state of the linear or of 
the zigzag configuration, as shown in Fig.[T|a). The pulse 
performs a quantum quench by bringing the central ion 
into a superposition of ground and excited states. In fact, 



at a time t after the first pulse the state takes the form: 
m)) = ^={\g)\cj> g {t))+^\e)\ct>S))), (2-1) 



where ip is a controllable phase and 

|^(t))=exp(-i£T a t/ft)|0(O)) 



(2.2) 



with s = g, e and H gi H e the Hamiltonians for the 
external degrees of freedom, accounting for the state- 
dependent potential. The free evolution is pictorially 
shown in Fig. [ljb) and leads to entanglement between 
internal and external degrees of freedom. After the sec- 
ond pulse, which performs a —tt/2 rotation of the dipole 
as sketched in Fig.[TJc), the probability of measuring the 
central ion in statefg) reads 

V g ^)= 1 -{l + Re[e^O(t)]}, (2.3) 



where 



0(t) = (Ut)\Mt)) (2.4) 



is the overlap between the two motional states. The con- 
trast of the Ramsey fringes is given by 



V(t) = \0(t)\ , 



(2.5) 



and depends on the time t elapsed between the two Ram- 
sey pulses. We note that the visibility of the interfer- 
ence signal is directly related to the Loschmidt echo, 
frequently used to describe the loss of coherence as a 
consequence of the interaction between a system and its 
environment |14) . 

Figure [2] displays the visibility of the interferometric 
signal, given by Eq. (2.5), as a function of the time 
elapsed between the two pulses. The visibility is eval- 
uated using the formula derived in Sec. |III[ The three 



plots correspond to three regimes we consider in this pa- 
per. In[2^a) the equilibrium configuration of the crystal is 
always a linear chain, which experiences a tighter poten- 
tial when the central ion is excited. Therefore the first 
quench does not change the equilibrium positions but 
rather the frequencies of the normal modes. The corre- 
sponding visibility exhibits sinusoidal oscillations and is 
close to unity. Figure [2^b) shows the case when the equi- 
librium configuration is a zigzag if the central ion is in 
the ground state, while it is a linear chain when the cen- 
tral ion is excited: The visibility decays quickly to zero, 
in agreement with the theoretical predictions for the de- 
coherence of a spin coupled to an environment close to 
criticality 4, 5 , but exhibits revivals with different peak 
heights. Figure [2]c), finally, corresponds to the situation 
when the equilibrium configuration of the crystal is al- 
ways a zigzag, which experiences a shallower potential 
when the central ion is excited. Here, the quench is also 
associated with a displacement of the equilibrium posi- 
tions. Similar to case (b), the signal decays and exhibits 
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(a) 



(b) 



(c) 




FIG. 2. (Color online) Visibility signal as a function of the 
time t elapsed between the Ramsey pulses for three 9 Be + 
ions with an axial trap frequency of v x = 2-zr x 1 MHz and 
transverse frequencies: (a) Vy = 27r x 1.565 MHz, (b) v y = 
2tt x 1.545 MHz, (c) Vy = 27r x 1.470 MHz. The critical value 
is v c = 2n x 1.549 MHz. The frequency of the transverse 

potential for the central ion when it is excited is ^Jvy 

where l>dip,g = 27r x 245 kHz. 



+ v. 



dip,!/' 



revivals. The decay of the visibility in (b) and (c) arises 
from the entanglement of the spin excitation with the 
crystal degrees of freedom, and the revivals are a signa- 
ture of quantum coherence that is stored in the whole 
system. Similar revivals have been experimentally ob- 
served for a single trapped ion [15j . They are here an 
intrinsic property of the many-body system and, as we 
will argue in the following, are a scalable feature that 
appears close to criticality. 

The details of the model that determine the proper- 
ties of the overlap integral, and thus of the visibility, are 
reported in the next Section. We note that in Eq. (2.4) 



we assumed that there is no mechanical effect associated 
with photon absorption and emission. This is the case 
when the internal transition is excited by means of a 
radio-frequency field [16] , or by a Raman transition with 
co-propagating beams [17] • The mechanical effects can 
also included in our formalism, see for instance Ref. [5], 
but will not change substantially the results for the cases 
illustrated in Figs.[2]jb) and (c). If the equilibrium struc- 
ture is a linear chain independently of the internal state 
of the ion, as in Fig. [2^a), a momentum kick would in- 
duce an oscillation about the equilibrium positions that 



would modify the signal. We refer the reader to Ref. [5], 
where a similar situation was studied. 



III. THEORETICAL MODEL 

We now give the detailed form of the Hamiltonian H 
that determines the evolution of the system. In the fol- 
lowing we will restrict the motion of the crystal to the x-y 
plane assuming a tight confinement in the z direction, so 
that the motion along this axis can be considered frozen 
out. The coordinates r = (x,y) give the position in the 
plane z = 0. This assumption is made for convenience: 
The calculations of this paper can be straightforwardly 
extended to three dimensions. 



A. Hamiltonian 

We first consider the internal degrees of freedom. We 
shall assume that only the central ion can be excited, 
while all other ions remain always in the ground state. 
The internal dynamics between the pulses can be re- 
stricted to the central ion, with Hamiltonian: 



H c \ = hhj. 



eg 



(3.1) 



where 



J eg 



is the transition frequency. The 



pulses are applied at time t = and t and correspond to 
a unitary operation given by the Pauli matrix <j x . 

The Hamiltonian for the external degrees of freedom 
of the ions, H mot , depends on the internal state of the 
central ion. We denote by the position and by p; the 
canonically conjugate momentum of the ion labelled by 
i. The corresponding energy reads 



-Hrnot = -ffkin + K 



'pot 



V c 



oul • 



(3.2) 



where -ffkin 



Sili Pi/(2 m ) is the total kinetic energy, 

N N 



Vcou\ — ^ 



1 



=i i=i 



4ire \n - rA 



(3.3) 



is the Coulomb repulsion, while the energy associated to 
the external potential takes the form 



Vr 



pot 



N 



i=l 



trap 



(ri) + Vdi p (r ie )|eXe| 



Here, the trap potential is 



nap(ri) = — {ylx\ + Vyyf) 



(3.4) 



(3.5) 



with v x , v y the trap frequencies along the axes x, y re- 
spectively, and the spatial part of the spin-dependent po- 
tential reads 



-^dip?/i c 



(3.6) 
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where the subscript i c labels the central ion. We as- 
sume i/dip is small compared to v y . The total Hamilto- 
nian which governs the dynamics between the laser pulses 
takes then the form 



H — H c \ + i7kin + Vt 



pot 



V Cc 



(3.7) 



In particular, the Hamiltonian H s = (s\H mot \s) deter- 
mines the dynamics of the external degrees of freedom 
when the central ion is in the internal state \s = g, e). 



The Hamiltonian of the crystal conditioned to whether 
the central ion is in state s = g,e takes the form 



H. 



is) 



off 



2A r 

V 



Pi 



^ 2m 



where V s is defined as 

a 2 

v fc = — 

3 drjdrk 



2N 

E 



in 



pot 



V c 



oul) 



(3.9) 



(3.10) 



B. Spin-dependent crystalline structures 



and V^ ot = (s\V P ot\s 



Hamiltonian (3.9) is transformed into a set of uncou- 



pled oscillators by an orthogonal matrix 1VF such that: 



We shall consider that the ions vibrate about their clas- 
sical equilibrium positions, with displacements from equi- 
librium that are much smaller than the inter-particle dis- 
tance p.8 . This situation can be achieved by laser cooling 
a hot cloud of ions confined in an electromagnetic trap, 
e.g. a Paul or Penning trap [IHlllOJ.The spin-dependent 
potential can be an optical potential, like a tightly fo- 
cussed laser beam propagating along and aligned with 
the chain axis, as discussed in Since this potential 
depends on the internal state, so does the crystal equilib- 
rium structure. For a fixed number of ions N the relevant 
parameters controlling the structure of the crystal are the 
aspect ratio a — v 2 /v 2 and the state-dependent shift to 
the aspect ratio a dlp — jv^,, where we consider that 
the spin-dependent force steepens the potential for the 
central ion. When all ions are in the ground state and 
a is larger than a critical value a c (N), the linear chain 
is stable, while at a c (N) it undergoes a continuous tran- 
sition to a zigzag [§[ [5TJ [22] . We shall assume that a is 
close to this critical value, so that the equilibrium struc- 
ture depends on the internal state. To study quenches 
across the phase transition by exciting the central ion, 
an accurate knowledge of the dynamical properties of the 
crystalline structures in both configurations is necessary. 



E 

jk 



m (wf ) 2 5 ln , 



where u)f are the normal modes frequencies and the cor- 
responding coordinates are related to the original dis- 
placements by the transformation Qf —J2k ^-kilk' with 
1 = 1,..., 2N. The second-quantized form of the Hamil- 
tonian is found introducing annihilation (creation) oper- 
ators b° (bf), with b s = yjmu?/(2h) [Qj + \Pj/(mu})\ 
and [b s j ,bf]=S jl : 



H 



0) 
cir 



2N 

EE 

s=g,e j = l 



s)(s\huj s . 



(3.11) 



The eigenstates are the number states {\n\, . . . , ^2Ar) s } 
with bfb s |ni,...,ri2iv) s = nj |m • • ■ti2n) s and Uj = 
0, 1,2,.. ., which form a complete and orthonormal basis 
for fixed s. The eigenstates of and are related 
by a transformation which is specified below and will be 
needed in order to study the dynamics of the system after 
the quench. 



D. Mapping between the normal modes of two 
different crystalline structures 



C. Spin-dependent normal modes 

In the following we introduce the notation for the 
normal modes of the state-dependent ion crystal, using 
s = g, e to indicate the internal state of the central ion. 
We denote by rf the equilibrium position of the i-th ion 
in the crystal for each internal state. A Taylor expan- 
sion of the potential V po t + Vcoui is performed to second 
order for small displacements qf around the equilibrium 
positions, qf = — r| . For convenience we will use the 
notation qi_ x t— > qj with j = 1, . . . , N and qi. y H ^ qj with 
j = N + 1, . . . , 2N. The following relation links the dis- 
placements between the crystal with the central ion in 
state g and e: 



q 9 = r e 



In order to evaluate the visibility, which is found from 
Eq. (2.4), we need to determine the transformation which 



connects the quantum states of the linear and zigzag 
structures. In this subsection we show that this is sim- 
ply found from the transformation which connects the 
ground states of the linear and of the zigzag configu- 
ration. This transformation is derived below, the final 



result is given in Eq. (3.35) 



For this purpose, we first consider the mapping relating 
the normal modes with displacement Q\ and Qf . This is 



found starting from Eq. (3.8) and rewriting it as 



(3.12) 



where d a - = r| 



r, is the difference between the equilib- 



(3.8) rium values for the coordinate . Inserting the definition 
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of the normal modes one finds 

Q 3 j =zZ T ^ + D h ( 3 - 13a ) 

k 

P?=J2?jkPk, (3.13b) 

k 

with Tji = J^k MjyM|j an orthogonal matrix and D 9 — 
J2k Mfej^fc the mode displacements. The transformation 
of the corresponding normal-mode annihilation and cre- 
ation operators is given by a Bogoliubov transformation, 
obtained by inserting the definitions of the operators into 
relations ( 3.13| , and which takes the form: 



k k 



(3.14a) 
(3.14b) 



ft fc 
Here, the real dimensionless coefficients Uj k , Vj k read: 



Ujk 



v jk 



T jfc 
2 

Tjfc 
2 



and fulfill the equations [53] 

£ ( U jk u lk - VjkVlk) 
k 



Vj,L 



(3.15a) 
(3.15b) 

(3.16a) 
(3.16b) 



Coefficient f3 9 describes a displacement of the corre- 
sponding normal mode: 



01 = 



— J - D 9 . 
2h J 



(3.17) 



After having obtained these relations we can now 
identify the transformation connecting the basis states 
{\m,n2, ■ ■ -) e } and {\ni,ri2, ■■■}„}■ Since every state 
of the bases can be generated from the corresponding 
ground state by applying repeatedly the correspond- 
ing creation operators, it is sufficient to find a map- 
ping between the ground states |0,0, . . . ,0) e = |0) e and 
1 0, 0, ... , 0) = |0) . Such mapping is given by a unitary 
transformation U such that 



|0> o =^|0) c 



(3.18) 



Operator U connects two Gaussian states and can thus 
be written as: 



U = 2?(7i, . . . ,72jv) S(fii,£L2,.. 



2N2N 



(3.19) 



where I? is a displacement operator and 7j are real 
scalars, while S is a squeezing operator that takes the 
form 



S = exp ^ " ^W^j > 



(3.20) 



with squeezing parameters £j k to be determined. With 
the help of the disentangling theorem, Eq. (3.20) can be 
recast into the convenient form 



C ry a b -A 1 ' 
o = Ze e e 



where Z is a scalar while 



et 
k ) 



jk 



jk 



(3.21) 

(3.22a) 
(3.22b) 



are operators, with Ajk a symmetric matrix. Details of 
the derivation are provided in Appendix [X] Application 
of operator (3.211 to the state |0) e gives 



|0) s = W|0) e , (3.23) 

with the non-unitary operator W defined as: 

W = ZV{ lll ... ll2N )e K . (3.24) 

We note that this operator was first introduced in [24] 
for evaluating the thermodynamics of interacting con- 
densates. 

We now determine the coefficients Ajk, the displace- 
ments 7j and the normalization constant Z . For this pur- 
pose we make use of relation bj \0) g = which must hold 



for any mode j of . Using Eqs. ( |3.14a| ) and (|3_23j, 

one obtains: 



= b 9 \0) g = 



^(u jk b%-v jk b%^+(3: 



L k 



W|0) e 



(3.25) 



which can be recast in the form 



l k 

E ( u jk - Vjk) ik + ft: 



|0) =0 (3.26) 



The latter equation has been derived from (3.251 multi- 
plying by 1 = WW -1 on the left side and making use of 
the relations 

W-% W = b% +J2 A kibf + lk , (3.27a) 



w- 1 bpw = bp + lk . 



(3.27b) 
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Equation ( 3.26[ ) is equivalent to: 

^ UjkA k i - Vji = 

k 

£(%* - v 3k )ik + 0° = , 



(3.28) 
(3.29) 



which must hold for all j,l = l,..., 2N. From Eq. ( |3.28[ ) 
one finds the coefficients 



a vik ■ 



(3.30) 



where one sees that A is real, with symmetry following 
from (3.16b), while from Eq. (3.29) obtains 



7j 



(3.31) 



where /3| has been defined. Finally, the constant Z is 
found from the condition that the norm of the ground 
state \0) g must be unity, g (0|0) g = 1, giving 



oo oo l+tl A m 



!>v t w|o) £ 



EE 

\n— m— 



AT A 

).! m! 



which leads to: 



Z = dot 



(l-A 



2\l/4 



|0) e , 
(3.32) 

(3.33) 



(details are given in Appendix |B[) . 
Using this result, we find: 

|0) =Z2? e (/3 e )e A |0) 



(3.34) 



with Ajk, Pj, and Z given in Eqs. (|3_30f, ( |3.3l[ ) 
and (3.33), respectively. Using this formalism we can 



now evaluate the visibility (which we defer for Sec. HIE) 
as well as the overlap between the two ground states: 

£o= e (0|0) g = Z e (0|P(/3 e )e A |0) e 

= Z exp | \ £ A^l }ex P {-l£ |^| 2 j . 

(3.35) 



Before we conclude this section, we also give the form of 
the squeezing parameters in operator S (3.21): 



£j* = E Aj* arctan h( a i)A* 



(3.36) 



where A-,-; is the orthogonal transformation diagonalizing 
A and a; are the corresponding eigenvalues. The param- 
eters £jk are real, since Ajk is real and symmetric. The 
derivation of Eq. (3.36) can be found in Appendix [A| 



E. Evaluation of the visibility 

We now derive an analytical expression for the visibil- 
ity. Our starting point is the overlap as a function of the 
time t between the pulses, as given in Eq. (2.4). Using 



- U s (t) |0) , we rewrite it as 



exp{H^t/ih} \0) g 



o(t) 



\UUt)Ue(t)\0) 9 



U e (t) \% ■ (3-37) 



We note that this expression is given up to a time- 
dependent phase, which depends on the difference be- 
tween the (classical) ground-state energies of the two 
equilibrium configurations |10j . Since this factor is ir- 
relevant for the visibility, it will be omitted from now 
on. Using Eq. (3.34), expression (3.37) can be cast in the 
form: 



0(t) = Z 2 



,.A f 



Vl(P e )U e (t)V e (ne A \0) c 



Z 2 e (f3 e \e Ai ^U e (t)e A ^\(3 e 



(3.38) 



where |/3 e ) e = T> e {j3 e ) |0) e and from the first to the second 
line we employed the relation 



V e ((3 e )e A = e A ^2? e (/3 e ) . 



with 



(3.39) 



(3.40) 



Ok 



Using the overcompleteness of the multimode co- 
herent states, the identity operator reads 1 = 
^sw / d 2 ai . . . d 2 «2Ar |a) e (a|, and Eq. (3.38) takes the 
form 



0(t) 



TT 



z*_ 

2N 



TT 



2N 



d ai . . . d 2 a2Ar 

(/3 e |e At ^C/ e (i)|a) e (a|e A ^)|/3 e ) ( 
r 

I d 2 ai . . . d 2 a 2 N 
Here, we have defined 

, Z^j \ 2 T 2 



Ca,f) 

and used 



-«*/?) 



(3.41) 



(3.42) 



e (f3 e \ e At ^ \a(t)) e = eVW-W^Fl a(t)) e , (3.43) 

where \aj(t)) = \otj exp{— iw|i}) is the time-evolved 

coherent state and /(a) = ^J2jk A jk a j a *k a 2./V- 
dimensional complex-valued function. The evaluation of 



Eq. (3.41) is just a matter of algebra and is shown in 
Appendix [Cj The result reads 



0(t) 



k-w 1 n _1 w 



Vdctft 



\Qo\ 



(3.44) 
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where the complex symmetric 4N x AN matrix fl and the 
4V-dimensional vector w are defined as 



CI = 
Here, 



and 



1 - A+ -iA~ 
-iA~ 1 + A+ 



A% = U jk [e-^^±l] 



-iS~ 



S j {F*)±S0 e )e-™P.. 



(3.45) 

(3.46a) 
(3.46b) 

(3.47) 



a linear (g > 0) or a zigzag chain (g < 0). The vertical 
axis gives the dimensionless parameter A, defined as 



A = 



'dip 

„2 



and related with the change in the potential on the cen- 
tral ion when it is in state \e). The equilibrium con- 
figuration when the ion is in state |e) is shown in the 
diagram as a function of g and A. We restrict to the 
case A > 0, consistently with the choice that the spin- 
dependent force is restoring. The ion in the excited state 
feels thus a change of the transverse potential correspond- 
ing to a vertical shift in the diagram as sketched by the 
green arrow. 



Equation (3.44 1 determines the visibility, that is plotted 
in Sec. IV for several parameter regimes in which this 
formula is valid. 

The short-time behaviour of the visibility is found by 



performing a Taylor expansion of (3.371, and reads: 
V(t) « 1 + f]t 2 /2 . 



(3.48) 



Quantity rj < 0, denoted in Sec. IV as the curvature, is 
determined by the variance of in the initial state: 



1 



(3.49) 



The functional dependence of r/ on the parameters in 
Eq. (3.44 1 is derived and given in Appendix [C| 



IV. QUANTUM QUENCHES AT THE 
LINEAR-ZIGZAG TRANSITION 

We shall now examine the visibility of the interfer- 
ometric signal when the chain is close to the linear- 
zigzag instability. We assume that a quench is per- 
formed by exciting the central ion in presence of a spin- 
dependent force. Due to the long-range interaction, the 
force on the central ion can induce a change of the equi- 
librium configuration of the entire crystal. In particular, 
if i/y < < Vy + v^ip y i the two equilibrium configura- 
tions corresponding to the different internal states can be 
a zigzag and a linear chain, respectively, provided that 
the correlation length is larger than the size of the sys- 
tem |10l I25j . The equilibrium configurations correspond- 
ing to the central ion excitation are represented in the 
diagram of Fig. [3] Here, the horizontal axis gives the 
dimensionless parameter g, which is defined as 



This parameter determines whether the equilibrium con- 
figuration corresponding to the central ion in state \g) is 
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FIG. 3. (Color online) Phase diagram for three ions as a func- 
tion of the dimensionless parameters g and A. The solid red 
line separates the parameter regions where the ions form a 
zigzag (bottom left) or a linear chain (top right). The crys- 
talline structure corresponding to state \g) is at A = 0. Three 
regimes are considered, depending on where the state before 
and after the pulse are located in the diagram: when they are 
both in the left region, g < g c , the chain is in a zigzag struc- 
ture; when they are both in the right region, g > 0, the chain 
is linear, while when they are across the line, g c < g < 0, the 
initial ground state is a zigzag and the final is a linear chain. 

Three situations will be discussed in the regime close 
to the linear-zigzag instability, indicated by the solid line 
of Fig. [3] The first one corresponds to the case in which 
the crystal is initially forming a linear chain (g > 0). 
The spin excitation then does not change the equilib- 
rium configuration, nevertheless it modifies the frequen- 
cies of the normal modes. An example of the visibility 
one measures in this case is displayed in Fig.^fa). When 
g < 0, the equilibrium configuration of the initial state 
is a zigzag. Whether the equilibrium configuration of 
the excited state is a zigzag or a linear, depends here on 
whether the shift A is below or above the instability line. 
For a fixed A, this defines a critical value <? C (A), such 
that at g — g c (A) < the crystal equilibrium configu- 
ration for the excited state is exactly on the instability 
line (note that g c depends on the number of ions N). If 
g < g c (A), hence, the crystal equilibrium configuration in 




FIG. 4. (Color online) Curvature rj of the visibility signal at 
short time, Eq. (3.481, as a function of g and A. Subplot (b) 
shows the corresponding contour plot. The curve has been 
evaluated for three 9 Be + ions with axial trap frequency v x = 
2tt x 1MHz. 



the excited state is also a zigzag (with however different 
transverse displacement as the initial one). An exam- 
ple for the visibility found in this case is shown in[2]Jc). 
If > g > <7 C (A), instead, the crystalline structures of 
ground and excited states are a zigzag and a linear chain, 
respectively. The quench hence drives the chain across 
the critical point, and a typical visibility signal is shown 
inFig.gb). 

The behaviour of the visibility for short times is char- 
acterized by a decay with quadratic dependence on the 
elapsed time, according to Eq. (3.481. Figure [4| a) dis- 
plays the parameter 77 as a function of g and A. Decay is 
faster in the region where the quench is performed across 
the phase transition, where the overlap between initial 
and final states is small. The plot is reminiscent of the 
features of the stability diagram in Fig.|3j as is visible by 
inspecting the contour plot in Fig. Hb). 

We now analyse the behaviour for long times. Figure[5] 
displays the density plot of the visibility as a function of 
the rescaled aspect ratio g and of the time t between 
the pulses. Three distinct behaviours are observed corre- 
sponding to the three regimes we identified. For g < g c 
the appearance of the revivals is periodic. The corre- 
sponding period diverges as g approaches g c . Each curve 
indicating a maximum of the visibility, moreover, ex- 
hibits an additional modulation, showing that the height 
of the revival peaks varies with g. The inset shows the 
behaviour at g — g c : Here, several peaks of the visibil- 
ity appear at short elapsed times, with rapidly vanishing 
height. In the interval g c < g < the periodic structure 
of the revivals is also observed, with decreasing period as 
g approaches 0. Finally, for g > the visibility is close 
to unity and exhibits some modulation for small positive 
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FIG. 5. (Color online) Density plot of the visibility signal 
as a function of g and t (in ^s) for A = 0.025. The other 
parameters are the same as in Fig. [4] The inset enlarges the 
region about g c for short elapsed times. 
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FIG. 6. (Color online) (a) Times (in fis) corresponding to 
the first peak of the visibility as a function of g for different 
values of A = 0.025 (symbol: A), 0.015 (□), 0.005 ( ). The 
dashed vertical lines denote the different values of g c - Inset: 
the same curves are plotted as a function of g — g c (A). The 
other parameters are the same as in Fig. El (b) Same as (a) 
but for a number of ions TV = 3 (O), 5 (0), 7 (A), 9 (□), 
11 (x) with A = 0.025. In the inset the curves are plotted as 
a function of g — g c {N) with the time in units of l/v c (N). 



g, with an amplitude that vanishes as g increases. 

Let us now make some considerations. In the first 
place, the value of g c depends on A. Nevertheless, the 
behaviour found in Fig. [5] is encountered for different val- 
ues of A, as is visible in Fig. [6j where we show the re- 
vival times at which the visibility is different from zero. 
This behaviour is also to large extent independent of the 
number of ions composing the crystal, as is indicated by 
Fig. [6j Here, one observes that all curves giving the first 
revival time for different numbers of ions exhibit a similar 
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functional dependence as g approaches g c (which depends 
also on the number of ions N) . This behaviour becomes 
evident by appropriately rcscaling the curves as shown in 
the inset. The width of the peak at g c in turn decreases 
as the number of ions is increased. 

Further information is gained by inspecting the Fourier 
transform of the visibility and of its logarithm, respec- 
tively defined as 

H"n) = ^£ dtV(t)e-^, (4.1) 

JiK) = ^ T dt ln[V(i)] e--"' , (4.2) 

where uj n — 2nn/T with n £ No and T is a time inter- 
val such that XVmin 3> 1, with hv m i n the smallest gap in 
Hamiltonian . Figure [7] displays the spectra corre- 
sponding to the visibility as a function of time in Fig. [5] 
and its logarithm. For g > 0, shown in (a), the main 
peak is located at twice the eigenfrequency of the zigzag 
mode, while for g < it is at the eigenfrequency of the 
corresponding lowest frequency mode, see (b) and (c), 
which becomes unstable when the critical value is ap- 
proached. The behaviour for g > hints to the presence 
of squeezing, originated by quenching the trap frequency 
(and thus the normal mode frequencies) by exciting the 
central ion. The frequency of oscillations of the visibility 
observed in Fig. [5Ja) corresponds indeed to 2u>f, which 
is twice the frequency of the zigzag mode of the linear 
chain. For g < 0, the main peaks are also associated 
with the vibrational mode that becomes unstable at the 
critical point, and which is the one that is most signif- 
icantly excited by the quench. The main peak is now 
at tuf instead of 2ujf because the dominant effect of the 
quench is the displacement of the equilibrium positions. 
In (b) and (c) minor peaks are present at the eigenfre- 
quencies of the modes which couple to the zigzag mode 
and in (b) also at sums of eigenfrequencies. 

The spectrum of the logarithmic signal as a function 
of g is displayed in Fig. [SJfor g < and A constant. One 
observes a main peak corresponding to the mode whose 
frequency goes to zero when g approaches g c and the 
crystal in the excited state becomes unstable. Close to g c 
one observes two signals which become more visible, that 
are at twice the zigzag eigenfrequency and at the sum of 
the zigzag mode frequency with the axial breathing mode 
frequency. They hint at the presence of single-mode and 
multimode squeezing due to the quenching, suggesting 
that some entanglement between the modes is generated 
by the quench. 

We now summarize our findings. In the first place, 
the visibility decays fast to zero when the quench is per- 
formed between motional states whose classical equilib- 
rium configurations differ. The decay is faster when the 
two structures have different symmetries, and thus differ- 
ent spectral properties, such as when the quench is across 
the phase transition linear-zigzag. Nevertheless, the visi- 
bility exhibits revivals as a function of the time after the 



quench. These revivals occur at the frequency of the low- 
est normal mode of the zigzag structure, which becomes 
unstable at the instability and corresponds to the zigzag 
mode of the linear chain. This mode is in fact the one 
which has the largest overlap with the difference between 
zigzag and linear structures. This feature of the visibility 
is scalable, the periodicity of the revival remains in fact 
finite as the size of the chain is increased. The height 
of the peaks, however, decays as N is scaled up, con- 
sistently with the fact that the amplitude of excitation 
of the zigzag mode by a displacement of the central ion 
decreases as the size of the chain is increased. In the 
thermodynamic limit, hence, the surviving feature is the 
decay of the visibility signal at short times, corresponding 
to the fact that the quantum superposition of the spin ir- 
reversibly dephases. When instead the quench is between 
two linear chains, there is a quasi-periodic rephasing of 
the system, corresponding to the creation of squeezing by 
sudden changing the trap frequency and thus the normal 
mode frequencies [55]. The amplitude of the oscillations 
decays to zero while the visibility reaches a constant value 
which approaches unity as g is increased. 

Before we conclude this subsection, we remark that the 
dynamics we consider in this paper is unitary: decay of 
the overlap signal is solely due to dephasing in the dy- 
namics of internal and external degrees of freedom, while 
external sources of decoherence and noise have been ne- 
glected. They are expected to introduce a damping factor 
in the overlap signal, setting an upper bound for partial 
revivals decreasing with time. Some estimates have been 
already provided in Ref. TT\ , where it was argued that for 
the parameters of existing set-ups the revivals of the vis- 
ibility should be observed. In particular, the main source 
of noise in ion traps is considered to be patch potentials 
at the electrodes [27J 121] • The measured heating rates 
depend strongly on the actual set-up of the experimen- 
tal apparatus, the smaller heating rates which have been 
reported correspond to timescales of the order of several 
milliseconds, which would allow one to observe several 
revivals of the visibility. The other important point re- 
gards the assumption that the chain is initially in the 
ground state of the vibrational excitations. Ground state 
cooling of ion chains composed to up to 4 ions have been 
successfully demonstrated in [29] , The visibility signals is 
degraded as the temperature T of the chain is increased. 
The functional dependence of the visibility on T is sub- 
ject of ongoing studies. 



V. CONCLUSIONS 

The dynamical properties of an ion crystal after a 
quench have been theoretically investigated, when the 
quench is performed by creating coherent superpositions 
of motional states close to and across the linear-zigzag 
structural transition. These dynamics have been related 
to the visibility of the signal when Ramsey interferom- 
etry is performed on one ion of the chain. The visibil- 
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FIG. 7. (Color online) The spectrum of the logarithm of the visibility, Ji(u„) in Eq. (4.2 1, corresponding to the plots in 
Fig. [2] The vertical straight red lines indicates the normal-mode frequencies of the crystal when the central ion is excited, the 
dashed green lines give the doubled frequencies, the da sh-do tted magenta lines show sums of frequencies. The insets show the 
corresponding spectrum of the visibility, J-(uj n ) in Eq. (4.1 1. The parameter g is, from left to right, 0.02, — 0.005, and —0.1, 
while the critical value is g c = —0.0165, and A = 0.025. 
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FIG. 8. (Color online) Density plot of the logarithmic spec- 
trum J-l(u) of the visibility for the parameters of Fig. [4] 
A = 0.025, and g < 0. The dashed line shows the critical 
value q c = —0.0165. 



ity decays at short times as the internal state becomes 
entangled with the motional state of the crystal, but ex- 
hibits periodic revivals at longer times, determined by the 
frequency of the zigzag mode. Further periodic signals 
appear at multiples of the zigzag mode and at sums of 
different motional excitations, suggesting squeezing and 
entanglement in the vibrational motion generated by the 
quench of the trap frequency. These spectral proper- 
ties persist as the number of ions increases, even though 
the heights of the revivals decrease. These results are 
based on a theoretical model which we report in detail 
and which allows one to calculate the visibility for dif- 
ferent parameter regimes. This model is valid as long as 
as the harmonic theory of the crystal is applicable and is 
thus reliable for the parameters we consider in this paper. 

Our analysis shows that, if the crystal is initially in the 
motional ground state, these features can be observed for 
parameters that are consistent with ongoing experimen- 
tal work. The signal is however degraded as the tem- 
perature is increased, the functional dependence of the 
visibility signal on the chain temperature is object of on- 
going studies. 

We conclude by observing that the visibility signal al- 
lows one to study the behaviour of the soft mode across 



the classical phase transition. Extensions of these stud- 
ies to the parameter regime where quantum effects at the 
phase transition are relevant [301 131] would allow one to 
extract the corresponding quantum fidelity, in the spirit 
of the work performed in |32j . and will be subject of fu- 
ture studies. 
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Appendix A: Multimode Squeezing Operator and 
disentanglement theorem 



In order to obtain Eq. (3.211 from Eq. (3.20) we fol- 
low the method of Bogoliubov and Shirkov [33 . The 
basic idea is best understood by first considering opera- 
tor e A( - cr++0 '-\ with A scalar, a z the Pauli matrix and a± 
the raising and lowering operators for a spin 1/2, such 
that [<7+,er_] = <J Z , [oz,e>±] = ±2cr±. The disentangling 
formula reads: 



e A(<T ++( T_) = g, 



r_|_ tanh(A) ^— a z ln[cosh(A)] ^cr_ tanh(A) 



(Al) 



and can be obtained using the procedure sketched in 
Ref. |34j . Assuming A to be a continuous parameter, 
one makes the ansatz e H<r++°-) = e /(A)CT +F(A)e 9(A)<T - 
with F(A) an operator such that F(0) — 1, while /(A) 
and g(X) are analytic functions of A with /(0) = and 
g(0) = 0. Under these assumptions F(X) can be written 



F(X) 



e - /(A)<r+ e A(<7 + +<7_) e -g(\)<J- 



(A2) 
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We take the derivative of F(X) and obtain a first- 
order differential equation which contains all opera- 
tors. The contributions from a + and er_ cancel out 
by choosing /(A) = g(X) = tanh(A), so that F(X) 
f x [) [— a z In (cosh A)], hence demonstrating Eq. ( Al ). 
This procedure can be generalized to show the equality 

e f(at 2 -a 2 ) = e l^LAat 2 e -(ata+i)ln[co S hA] e -i^a 2 

where a, a T are the annihilation and creation operators 
of a harmonic oscillator. 

Moreover, we can use the procedure sketched above in 
order to disentangle the multimode squeezing operator: 



cxp j \ E fe^W - Zjk a j a k) | 



with 



% e k Ejfc A ik°-i W e ~J2jk B jk a/a k g - | C jk a^a k 

Ajk = y^tajih(x;)Aj;A fc j , 
I 

I 

C jk = J^tanhCxOAJjA;, = , 



Z =exp|-^iln(coshxi)| 



(A4) 

(A5) 
(A6) 
(A7) 

(A8) 



This can be done after observing that, since £ is com- 
plex symmetric, we can perform Takagi's factoriza- 
tion [35] £ = AxA T , where A is unitary and x — 
diag({xi, X2i ■ ■ ■ }) is diagonal with Xj > rea l an d 
non-negative (this factorization exists for any complex 
symmetric matrix). This defines the transformation 

bjt = J2jk A fej a fc t : ^ = J2jkKj a k f° r a new set of 
operators for which the squeezing operator is in diag- 



onal form, exp 



t2 



These new op- 



erators have bosonic commutation relations, [bj,bf.] — 
[bji,bk*] = 0, and [6^-, & fc + ] = J2i A ife A y = S jk since A is 
unitary. Therefore, operators of different modes factorize 
as EL ex P {if (bj^ 2 ~ bj 2 )}> an( i ° ne finally obtains 



GXP { 2 E fefc a / a fc f _ Cjk a 3 a k) 

= yj e ^»^^t2 e _ (i) .t 6 . + i )ln ( coshx .) 



e 2 j 



The terms belonging to different modes commute now, so 
bringing factors with operators 6j to the left and factors 
with b 2 to the right and writing them as a function of 



operators a k and a, , one obtains Eq. ( A4 1 . 



Appendix B: Calculation of the Normalization 
Constant 



In order to calculate the constant Z we use the nor- 



malization condition of the states as stated in Eq. (3.32 1, 

AT A 



' oo oo 

EE 



J ! TO ! 



|0) 



Since A contains only creation operators, only the sum- 
mands with m — n give a contribution: 



z- 2 = 



f oo » 4.n . n 

A T A 

^ 7^)2" 
\n=0 v ' 
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^ (n!) 

n=0 v ; 

where W n is defined as 
1 



E7^<°i At " A >> = XX< 



(Bl) 



n=0 



w n = 



{2 n nl) 2 



Ajl]2 ' ' ' A^n-ljan^fcl' 



kl---k2n 



■■■A 



• &L„ |0> - (B2) 



The sum contains (2n)! summands (which contain 2n 
Kronecker-delta symbols) which do not vanish, corre- 
sponding to the number of all pairs of sets of indices 
{ji ' ' ' 32n} a nd {k\ • • • k2 n } which are identical (apart for 
a permutation within the same set). For example, 

w i = p E A hh A kxk2{ S jxki^j2k2 + S jik 2 8j 2 kiY 

k±k 2 

while W2 has already 24 summands, we write only two 
of them exemplarily: 

y t Aj 1 j 2 Aj 3 j i Ak 1 k 2 Ak 3 k4 ' ^jik4,5j 2 k 1 5j 3 k 2 5j i k 1 



kxk 2 k 3 k t 



(B3) 



y ] Aj 1 j 2 Aj 3 j 4 Ak 1 k 2 Ak 3 k 4 ' ^jik 1 Sj 2 k 2 Sj 3 k4^j i k 3 



0l32] 3 3i 
k!k 2 k 3 k 4 



(B4) 



We now associate with each summand in W n a graph, 
which we call a n-graph. For this, let the first n coeffi- 
cients be represented by n pairs of adjacent circles in an 
upper row, while the second n coefficients are represented 
by the same number of pairs of circles in the lower row. 
The indices {j\ ■ ■ ■ j2n} s-nd {k\ ■ ■ ■ k2 n } are filled in in 
correct ordering into the circles such that there are only 
j's in the upper and only fc's in the lower row. 



j: 0© ®@ 



*: 0© 0® 
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Then for each Kronecker-J we need to connect the corre- 
sponding two circles by a line. We find easily that each 
circle must be connected with another, and that there is 
a total of 2n lines. Thus each circle has exactly one line. 
For instance, the graphs for (B3) and (B4) are given by: 




and 




respectively. If wc evaluate (B3), we find that it yields 



Tr(A 4 ), while JbS can be factorized into two terms 



32^2 



3132 
k 1 k 2 



3i k 3 



3334 
k' 3 ^'4 



which give [Tr(A 2 )] 2 . This factorization can also be 
shown graphically, 






0® 
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Thus, a graph may be decomposed into a product of 
fully connected subgraphs or clusters. An n-graph can 
be decomposed into a product of mi 1-clusters, m2 2- 
clusters, . . . , and m„ n-clusters, where the to; fulfill 



i=i 



mi I = n. 



(B5) 



The evaluation of an Z-cluster always yields Tr(A 2i ), and 
there are 2 l l\ 2 l ~~ 1 (l — 1)1 ways to draw an Z-cluster. So 
we are motivated to define the l-cluster integral by the 
sum of all possible clusters for I pairs of circles in each 
row, which after evaluation is given by: 



k = 2 l l\ 2 l - 1 {l - l)!Tr(A 2 



(B6) 



We have bo — 1 and b\ = 2 Tr(A 2 ), which finds its graph- 
ical representation by 



oo 
oo 



Accordingly one can draw the cluster-integrals for the 
higher orders. Here we have not filled out the circles, 
since the cluster integral is independent of the indices 
which are assigned to it. It is clear that for a given set of 
indices {i i , %i , . . . } the circles have to be filled in the same 
ordering in each summand, and without loss of generality 



one can fill the circles in the natural ordering ,j i2 ,...) 
and (fe^, ki 2 , . . .) where i\ < i% < ■ ■■ . The total set of 
indices cannot be split arbitrarily in between the clusters, 
since pairs of indices of the form (J21— 1,321) always belong 
to the same cluster. 

To resume the calculation, wc note that 



W n = 



1 



(2™n!) 2 



(B7) 



where S{mi} is the sum over all possible graphs described 
by the set of integers {mi}, and the primed sum denotes 
a restricted summation over all sets {mi} which fulfill 
equation (B5|. We see that 



s{m«} = x>r&: 



m 2 
2 



(B8) 



where the summation over V extends over all 
possible ways of distributing the two times n 
pairs of indices {(Ji'Jt), ■ ■ ■ , (jin-ilhn)} and 
{(kx; k-i), ■ ■ ■ , (fen-i! &2n)} i n to the circles obtaining 
only distinct graphs. So there are n! ways of distributing 
these pairs (the ordering of a pair is already contained 
inside the cluster integral). A permutation of two 
Z-clusters with the same I does not give a new graph, 
therefore we get a factor J^[ z (m;!) . Moreover, a per- 
mutation of pairs inside a cluster integral does not give 
a new graph either. Thus we get a factor IIz^O 2 ™'- 



Equation ( B8 1 is then given by 



S{mi} = (n!) 2 Y[ 



1 nil 



1=1 



mi\(l 



! (mam, 



(B9) 



Replacing in (B7) one gets: 



1 1 h m > 

V 1 {mi} '- 1 V ' 



1=1 

h 



— T'T\- 

2 2n ^ lL mi. 

{mi} £=l vv 

y' 1 frA 

A^i <Y2(mil+m 2 2+---) LL m \ 
{mi} 1=1 
n ^ 

e n 



hi 

w 



{mi} 1=1 



(2H 



(BIO) 



We can now insert this result in Eq. (Bl) and obtain 

k 



ee n m i 

n=0 m, 1=1 



{2H\f 



(BH) 



Summing over all {m{} followed by summation over all 
n is equivalent to summing over all mi, ni2, . . . from to 
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oo separately, so we can replace the restricted sum: 



oo oo oo ^ 



mi=0m 2 =0 8=1 vv y 

6/ ^ m '" 



n 

1=1 

oo 

J|exp 



y - , 

mi=0 1 vv 7 



8=1 



(2'/ 



H2 



exp 



E 

,8=1 



(2'Z!) 2 



(B12) 



Using equation (B6), we finally get 



exp 



E 

.8=1 



Tr(A 2 ') 
21 



exp 



Tr 



In 



1 - A 2 



(B13) 



which is valid if 1 — A 2 is non-singular. To show that 
this is true it is sufficient to show that any matrix norm 
of A is smaller than one. Using the spectral norm ||-||, 
the form of Eq. ( A7 1 , and the submultiplicativity of the 
matrix norm, we have ||A|| 2 < ||A|| 2 ||tanhx||2||A T ||2 = 
Ijtanhxlb- Using the fact that x 1S diagonal, real and 
positive, the spectral norm is equal to tangent hyperboli- 
cus of the largest eigenvalue of \- Thus ||A|| 2 < 1 as the 
tangent hyperbolicus is smaller than 1 in its full domain. 
Equation (B13l can thus be cast in the compact form: 



exp 



1 



Tr 



In 



1 

1 - A 2 



det 



(1 



i 2\ 4 



(B14) 



Appendix C: Calculation of the Overlap Integral 



We consider the integral (3.411 and first remove the 



time-dependent phase factors e 1Ujt from the integration 
variables otj in a*(a(t) — f3) by shifting it to the coef- 
ficients A jk by defining A jk (t) = A jk e' i{ui+Uk)t ■ We 
merge all terms into a single exponential whose exponent 
reads 



1 E ^ aj ^ ^ Ajh ^ ~ Sjk ^ ^ ak 



2 ^ \a 

jk 



Sjk A jk I \a k 



J2 Sj [/?>,- ~ E S i W'*"'**, + G* (p) + G(p) 



with 



and 



G(/3) = E^r^--E 



\pj 



2 j z— / 2 

jk j 



(CI) 



(C2) 



We now express the integration variables by their real 
and imaginary parts, ay = Uj + ivj and a* — Uj — iVj. 
The quadratic term is written as 



s jk - A+ - iA, fe \ / Uk 



S (7, 



-iA jk S jk + At k J V.'V. 



with complex symmetric matrices 



1 



a£ = 2(4**(*)±4**(°)) 



(C3) 



The linear term in the exponent takes the form 
-Ej [ s ? u j- iS I v i] with 



Sf = S j [/3*]±S j [l3]e- i ^ t . 



(C4) 



Introducing the vector w = (u, v) T where u — 
(ux, . . . , u 2 n), v — (wij • • • i v 2n), we can write the overlap 

as 



with 



-iS- 



1 - A+ -iA" 
-iA- 1 + A+ 



(C5) 



(C6) 



The result of the integral in Eq. ( C5 ) is y 
and Eq. (3.41) can be cast in the form 



<D(t) 



p 2Re{Go} e iw T n- 1 w 



with 



y/dctU 



Ajk ^p* \Pj I 



(C7) 



G =E^r^r-E 

jk j 



Using Eq.(3.35) in Eq. (C7), the visibility can then be 
cast in the form of Eq. (3.44). 



The convergence of the integral in Eq. (C5) is verified 
by showing that the matrix Q, = 1 — B, with 



B = 



A^ 
iA" 



iA" 
-A- 



(C8) 



has only eigenvalues whose real parts are greater than 
zero. For this purpose we consider the spectral radius 
of B, p(B) = max{|A B |}, where Ab is an eigenvalue of B 
and which fulfills p(B) < ||B|| for any matrix norm |35j . 
If 1 1 B 1 1 < 1 , it follows that all eigenvalues of B lie within 
a circle with radius p(B) < 1 centered around 1. Then, 
all the real parts of all eigenvalues of O = (l — B) are 
greater than zero. B can be brought to block-diagonal 
form Db by a similarity transformation with an orthog- 
onal matrix Mb: 



B = 



1 fl 



(k{t) 
I A 



(C9) 
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Thus ||B|| = ||M B D B Mg|| < ||Mb||||Db||||M b || by the 
submultiplicativity of the matrix norm. The spectral 
norm of the orthogonal matrices is 1, and the spectral 
norm of Db, ||Db || 2 = max {||A(i)|| 2 , ||A|| 2 }, but since 
||A(t)|| a - ||A|| 2 , wc have ||B|| 2 = ||A|| 2 . 
We now proceed to perform a Taylor expansion of 



Eq. (3.44) for short times. For this purpose we first bring 



expression (3.44) into a more convenient form, using the 
definitions 



3 = 1 + A+ . 



T = 1 - A+ 



(CIO) 



The determinant and the inverse matrix can be calcu- 
lated with the help of the corresponding identities for a 
partitioned matrix [36 , 



det Q, = det 3 • det 6 . 



(Cll) 



and 

o- 1 



e- 1 ie-^A-s- 1 

iS^A-G- 1 H- 1 - S- 1 A-0- 1 A-S- 1 ' ' 



where is the Schur complement of S given by 

9 = t + a-s^a - . 



(C12) 



The equations hold provided that 2 and are non- 
singular, which is true as shown in Appendix \K\ 
Expanding the overlap around t = 0, we find 



0{t) w 1 - iO x t - -0 2 t 2 . 



(C13) 



with O x = O(0) and 2 = O(0), which leads to the 
expression of the visibility in Eq. (3.48), where v = 

-(o 2 -ol). 
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